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pH I Abstract. We present a review of our recent works directed towards discovery of a periodic, 

^ ■ kink-like and soliton- like travelling wave solutions within the models of transport phenomena 

and the mathematical biology. Analytical description of these wave patterns is carried out 
I by means of our modification of the direct algebraic balance method. In the case when the 

analytical description fails, we propose to approximate invariant travelling wave solutions 
by means of an infinite series of exponential functions. The effectiveness of the method of 
approximation is demonstrated on a hyperbolic modification of Burgers equation. 
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^! 1 Introduction 

. It is well known that methods of obtaining general solutions for the majority of nonlinear evo- 

lutionary PDEs actually do not exist. The already known methods of integration, such as the 
inverse scattering technique, are applied merely to so called completely integrable equations 
These equations form a very important but relatively small family, which does not include 
e.g. nonlinear dissipative models, describing the spatio-temporal patterns formation and evolu- 
' tion [2]. One of the few alternatives to numerical methods of the investigations of those models 

that are not completely integrable are the symmetry methods which are effective when the 
system under consideration admits a non-trivial group of transformations. Yet, these methods 
cannot be treated as universal tools for obtaining analytical solutions. In fact, the self-similar 
reduction, when applied to a nonlinear PDE, gives another nonlinear differential equation with 
a fewer number of independent variables and it is difficult to expect that the procedure will 
lead to quadratures in every particular case. Obviously the chances to obtain this way the 
analytical description of solutions with the required properties and asymptotic behavior (such 
as e.g. periodic, kink-like or soliton-like travelling wave solutions) are, generally speaking, very 
small. 

In the case when the self-similar reduction leads to a system of autonomous ODEs the answer 
to the question on whether the solutions with the given properties do exist can be obtained by 
means of qualitative theory methods lEj , which allow for studying of whole families of solutions 
with the given symmetry. The qualitative theory methods are very effective if one wants to know 
whether the family of self-similar solutions possesses solutions with the given properties. Yet in 
many cases mere information about their existence is not sufficient. Very often one needs to 
have an analytical expression for the given solution and this is the case when it is necessary to 
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investigate the influence of small addends upon the stability and the evolution properties of the 
invariant solutions. 

Compared to classical symmetry methods, more powerful tools for finding out solutions with 
the given properties, can be provided by so-called direct algebraic balance method, based upon 
representing solutions in the form of rational algebraic combinations of properly chosen functions, 
closed with respect to the algebraic operations and differentiation [3IH1I11E1- However, within 
this method one must solve, as a rule, very large systems of nonlinear algebraic equations. 
Regardless of this purely technical reason, the ansatz-based method cannot be treated as the 
remedy, because there are known evolution PDEs which possess analytical solutions that cannot 
be obtained, at least, within a particular set of elementary functions (hyperbolic, trigonometric, 
polynomial, etc.). And it is rather doubtful that a fully universal ansatz for finding out the 
analytical solutions to nonlinear PDEs can ever be found. Nevertheless, it is possible to answer 
the question of whether the solution of this or that sort does exist within the given set of 
invariant solutions, whenever the procedure of a self-similar reduction leads to the system of 
autonomous ODEs. And in the situation when the direct algebraic balance method, based 
on the representing solution we are looking for in the form of finite combination of the given 
functions fails, we propose to express it by means of the infinite series. 

The aim of this study is to review our recent work directed towards the description of the 
already mentioned wave patterns, i.e. the kink-like and soliton-like travelling wave solutions, 
in various non-local models, simulating the transport phenomena. In our earlier work we in- 
vestigated merely the set of travelling wave (TW) solutions, being invariant with respect to 
the group of the spatio-temporal translations, nevertheless, systems obtained via the symmetry 
reduction prove to be sufficiently difficult to analyze and solve exactly. Therefore we use a 
complex approach including a qualitative analysis, our own interpretation of the ansatz-based 
method and approximation of the invariant solutions by means of infinite series of exponential 
functions. 

The structure of the article is as follows. In Section 2 we present a certain modification of 
the direct algebraic method, put forward in papers OIHlin]- Employing this modification we 
obtain the kink-like and soliton-like TW solutions for the hyperbolic modification of Burgers 
equation, telegraph equation, nonlinear d'Alembert equation and the somitogenesis model. In 
Section 3 we present the results of a qualitative investigations of a set of TW solutions for the 
hyperbolic modification of Burgers equation. It demonstrates the existence of the TW solution- 
like solutions, yet these solutions cannot be described by the finite algebraic combination of 
elementary functions. Therefore we propose to approximate the invariant solitary wave solutions 
by an infinite series of the exponential functions. 

2 Exact soliton-like and kink-like solutions 
for transport equations 

2.1 Modification of the direct algebraic method and its apphcations 

The essence of the unified direct algebraic method in E. Fan's paper [7], is based on the observation 
that particular solutions of any system of PDE's 



which do not depend on (t, x) coordinates in explicit form, can be presented as a polynomial 
expansion 




i = 1,2, ... ,m 



V = 1,2,...,?! 




s 




^ = X ^ vt 



(2) 



/i=0 
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where are unknown parameters, while function (/)(^) satisfies the equation 

m 



A 



v=0 



Note that in order to obtain non-trivial solutions, one should "balance" the numbers r and s 
(for details see 0). 

Depending on the conditions posed on the parameters Cy, solutions to the equation (jS)) are 
expressed by elliptic Jacobi functions, and in particular cases by hyperbolic or trigonometric 
functions [Jj. The properties of these functions are inherited by the solution of the initial 
system which can be presented in the form 0. In fact this methodology is constructive and 
algorithmic if the functions Hy (•,-,...) arising in are algebraic or rational ones. With this 
assumption, for any set of functions {it*}™^^^ that can be represented as @, the LHS of the 
equations can be presented in the form of finite linear combinations of functions (t>^{C) and 
4>\0m, A:,/ = 0,1,2,...: 

H,{u\ui,ui,...) = Y, buano + m E ^-^^^(0, = i, . . . , n. 
0=0 7=0 

When equating to zero all the coefficients of these decompositions, one obtains a system of non- 
linear algebraic equations, that determine the particular solutions of the initial system of PDE's 
(for more details see jjj). It must be stressed that the method does not guarantee obtaining 
of the general solution, even within the set (2), for its effectiveness depends on one's ability to 
solve a system of nonlinear algebraic equations. On the other hand, it enables calculation of 
particular solutions possessing some predictable properties. 

An enforced version of Fan's method was put forward recently in |S1 IH]. It was used for 
determining of exact solutions of the equation 

Ut + AUU^ - KUxx = f{u)- 

More precisely, the following ansatz was proposed: 



u 



40. 



^ = x + vt + xo, k = l,2,..., (4) 



k 

where z(^) = ^ 0'n(f>^{Cji 0(^) is a function satisfying equation Q. Due to this combi- 

nation, the multi-parameter families of the exact solutions were obtained in the situation when 
the pure form of Fan's methodology does not work IS]. 

On analyzing different versions of the ansatz-based method, one can conclude that effective- 
ness of their employment is based on the mere fact that the family of functions <j)^{0 and 4>'^4>{£,) 
is closed with respect to the algebraic operations and differentiation. In view of this, a quite 
natural generalization of the already mentioned ansatzes can be as follows: 

mi 1712 

,,,, Ea^.^'-io+mEbuno 
^^^^ EcxcpHo+mEd.no 

A=0 K=0 

where the function (p{^,) still satisfies the equation Q, but, in contrast to (jlj), dependence 
between the functions / and g is not assumed from the very beginning. Effectiveness of the 
ansatz (jSj) is demonstrated by applying it to some transport equations given below. 
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2.2 Soliton-like and kink-like TW solutions 
for the generalized transport equations 

Let us consider the following equations: 

TUtt + AuUx + But - KUxx =/("") = ^ XuU'^, (6) 

where r, A, B, k are non-negative constants. For ^ = equation coincides with the nonlinear 
telegraph equation; for j4 7^ it coincides with the hyperbolic generalization of Burgers equation, 
while for A = B = Q - with the nonlinear d'Alembert equation. The hyperbolic modifications 
of nonlinear transport equations arise naturally when the effects of temporal non-locality are 
taken into account jllj . 

Assuming that the solitons and kinks can be expressed by powers of function sech(,^), which is 
the particular solution of equation © and function sinh(^) which appears in the odd derivatives 
of the function sech(^), we use the following ansatz: 

mi m2 

a^sech^(a^) -l-sinh(a^) ^ busech^ [a^) 

"«) = il = ^? ? m 

£ c^sech^«) + sinh(QO £ d^sedf {a^) 

7=0 (7=0 

or, equivalently, 

m 

a/,exp(^aO 

n(0 = . (8) 

Inserting the ansatz © ( or (jH))) into © and executing all the necessary operations, we 
obtain an algebraic equation containing, respectively, functions sech^(a^), sech'^(a^) sinh(a^) 
or exp [/ia^] . Deeming them functionally independent and equating to zero the corresponding 
coefficients, we go to the system of algebraic equations. We do not demonstrate the details 
of these calculations since they are simple but cumbersome. The calculations were performed 
through the aid of a software application "Mathematica" . The results obtained are presented 
below. 

I. For arbitrary A, B and f{u) = Aq + Ain(t, x) + X2u{t, x)"^ + X-suit, xY the function 

u{t,x) = - ; r— TT (9) 

satisfies equation © if the following conditions hold: 



A 



-oooia + hQ)^ Ai = — ^ \abobi{BveA + /iG^) + AgooaiA^ 







^2 = T:T:^[ablbl{BvA + he) + XsA^e], a = [-2hablbl + X3A'] . (10) 
oqOiA^ aboOiA 

Here and henceforth we use the notation h = a{v'^T — k), A = aibo — aobi, @ = ai6o + 
oobi- Equation @ defines a kink-like regime when 69^2 > and ao/^o 7^ 02/^2- Using the 
conditions (|lUj) . we can express the unknown parameters from formula Q by the parameters 
characterizing equation yet, in the general case it is too cumbersome. It is much easier to 
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do when a = 2^J —\q\2/v^ oq = — ai = A0/A2, 5o = 61 = 1. From these conditions we obtain 
the solution: 



u 



where 



(t, x) = ^J-\q/\2 tanh 



a/— A0A2 



(x + vt) 



X2{AB + ^A^B'^ - 8KX3 + IGkA^t) A0A3 

, Ai 



2A3 - 4A^r 



As 



For K = l,r = 0,yl = 0,-B = l function Q coincides with the solution obtained in jH]. 
Another example of the kink- like solution defined by the formulae Q, ()lUjl is as follows: 



u{t, x) 



= X + vt, 



6o[l + exp(2aO]' 

where 60 = ( - A2 ± y^A^ - 4AiA3)/Ai, ^ = Aq = 0, a = -(2A2 + 3boXi) /[4Bvbo], and 

^ ±V;^(2A2 + 360A1) 

V4Air + 4boX2{B^ + 3Ait) + 6^Ai(2S2 + 9Air) ' 

II. For A = 0, arbitrary B and /(n) = Aq + A]^/2^(*5 x)2 + Aiu(t, x) + A3/2^^(i, x) 2 + X2u{t, x)^ 
function 



x) 



ao + aie' 



a{x+vt) 



bo + 6ie°(^+'^*) 

satisfies the equation (jH)) providing that the following conditions hold: 
Ao = 2alalah/A'^, Ai/g = -2ooaia(3/ie + BvA)/A'^, 



'1/ 

Ai = 2a(/i(3G2 - A^) + BvAQ)/A^, 
A2 = 66g6?/ia/A2. 



A 



3/2 



-26o6ia(5/iG + 5f A)/A^ 



This solution defines the solitary wave regime if > 0, |ao|/|6o| = 1 02 1 /I ^2 1 while ao/^o 7^ 
a2/fe2- 

III. For B = 0, arbitrary A and f{u) = Xiu{t, x) + X3u{t, x) function 



u{t, x) 



-of - 3af 026^ + 3aia^e2< + a^e^"? ' 



^ = X + fi, 



(11) 



satisfies (jSJ, when Ai, A3 are positive and the parameters are as follows: 60102 = — yAs/Ai, 
a = "v/AiAs/A, V = iby^]4^/A3 + K)7r. This solution is always singular, because for arbitrary 
values of the parameters the expression in the denominator of the formula Ullj) vanishes for some 

e e M. 

IV. Now let us consider the case A = B = 0. 

IVa. For f{u) = Aq + Xiu{t, x) + X2u{t, x)^ + A3n(t, x)'^ function 



x) 



ao 



bo + 26ie°(^+''*) + 6oe2°(^+''*) ' 
satisfies equation (jUJ when the following conditions hold: 

Ao = 00(2^0^0 — o-ibo — aoai6i)a/i/A^, 



(12) 
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= {{albl + 4aoai6o6i + al{-6bl + 6? ))a/i/A^ 
X2 = 3bo{2aobl-aiboh-aobl)ah/A'^, X3 = -2bo{bl - bl)ah/ A"^ . (13) 

For oo / 0, 60 / 0, |ai| + 7^ equation ((T^ defines the soliton-like solution. One of 
the parameters contained in (|12jl can be chosen arbitrarily, while the rest can be expressed, 
using ()13() . as the functions of the parameters, defining equation 0. We omit doing this in the 
general case, but present one particular example. Thus, for 61 = 0, 60 = 1 and arbitrary a we 
have the solution 

with 

ao = -A2/(3A3), ai = ^j2{Xl - X,\s)/Xl v = ±^[Xi - Xl/iSXa) + na^] /(ra^), 

and the additional condition Aq + AiAq + X2aQ + X^Gq = 0. 

VIb. For f{u) = Xi/2u{t, xfl'^ + Aiii(t, x) + A3/2u(t, xf^"^ + A2n(t, xf we obtain a soliton-like 
solution 

(ea(^+^'t) + 1)4 



with 

= -3a/i/^ii Ai = (1260 + 46i)q/i/^i, 
A3/2 = -(156g + 106obi)aV^i, A2 = (66^61 + 66i])a/i/^i- 

1 3 

Vic. For f(u) = Ao + X\i2u{t^x')'^ + Ain(t, x) + A3/2^i(i, a^) ^ a localized wave packet 
_ (aoe^"(^+^*) + (2ao + 4ai)e"(^+^^) + a^f 

defines a solution of if the following conditions hold: 

Ao = 2aQ(ao + ai)a/i/ai, A1/2 = (9^0 + 6aoai)a/i/ai, 
Ai = (12ao + 4ai)a/i/ai, A3/2 = —5ah/ai. 

VId. For /(n) = Ain(t, x) + X^/2u{t, x) ^ + A2?i(t,x)^ the function 
with 

Ai = Aah, A3/2 = — lOaia/i, A2 = (6af — 609)0/1 

defines a soliton-like solution of the equation Q. 
Vie. Finally, let us consider the equation 

TUtt - KUxx = Ao + AiU + X2U^ + XsU^. (14) 

Substituting the ansatz u = ^ = x + vt into the equation ((Tl)) . we obtain, after one 

integration, the following ODE: 

^ = ± V Co + ClU + 02^2 + CsU^ + C4U^, 

where cq is an arbitrary constant, ci = 2Xq/5^ C2 = Ai/5, C3 = 2X2 /{2>5), C4 = X3/{26), 6 = 
Tv^ — K. For this equation the classification given in [7j can be applied: 
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Figure 1. Phase portrait of dynamic system describing tire TW solution of system case IV. 

(a) if A2 = Ao = 0, then equation (|14|) possesses a soliton-like solution 

u{t,x) = ^-2Xi/\3sech{^/xij6C), Xi > 0, A3 < 0; 

(b) if Ao = A2 = 0, then equation (|14|) possesses a kink-like solution 

u{t, x) = V-A1/A3 tanh [^/^/(2(5)^] , Ai < 0, A3 > 0; 

(c) if Ao = A3 = 0, then equation (jTl]) possesses a soliton-like solution 

u{t, x) = - [3Ai/(2A2)] sech^ [^/xj6C/2] , Ai > 0. 

More general results can be obtained using the fact that H14|) corresponds to a Hhamiltonian 
case. Using the standard procedure, we obtain the following integral: 

C= / , 15 
J ^2H + 2(5-i(Aon + AiuV2 + Mw' + \j.u^l\) 

where H is the value of the Hamiltonian function on the given phase trajectory. The integral H15|) 
can be calculated explicitly in terms of elliptic functions for different values of H and as special 
cases contains all the solutions obtained in part IV. An example of phase portrait for A = 1, 
Ao = 2, Ai = —11, A2 = 15, and A3 = —6 is presented in Fig. ^ 

2.3 A model of mathematical biology 

Now we consider the following system |12j : 

ut + Auu^ - Uxx = fiu, z), 

zt = g{u,z), (16) 

where /(n, z), g{u, z) are given polynomials, and A is a positive constant. 

Due to the fact that the system is autonomous we apply variables in a travelling wave form: 

u{t,x) = V{0, z{t,x) = U{0, 

where ^ = x -|- /it, and we reduce system ()16|) to the following system of ODE: 

txV + AVV' -V" = f{U,V), ^iU' = g{U,V). 
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To simplify the computations, we restrict the family of admissible polynomials in RHS as 
follows: 

f(U, V) = aiU + a2V + a^UV, g{U, V) = hU + + b^UV. 
In order to obtain soliton- and kink-like solutions we proposed the following ansatz 



that shows the proper asymptotic behavior. 

The following kink-like solutions were obtained: 

f (?) = viO 



qi + q2e^ ' Aqi + Aq2e^ ' 

where 

ai=Ai-l, 02 = 0, 03 = , 61 = , 62 = 0, 63 = —. 

Pi 2gi 2 

For other family of parameters 

2 1 -2qin 2qi 

ai = n -1, ^2 = — i , 03 = , 

Api pi 

bi = , ^2 = /u , fea = -An 

2qi 

the ansatz provides the following soliton-like solutions: 
4pii?2e^ 



{qi+2q2e^r 

8qiq2e^ 



A{qi + 2(?2e«)(2g2(^ - l)e« + - 1)) ' 



3 Approximated soliton-like solutions 
of the generalized Burgers equation 

3.1 Existence of the soliton-like regimes within the set of TW solutions 

We showed in the previous section that the ansatz (jH]) is quite effective in obtaining travel- 
ling wave solutions with given properties and asymptotics. Yet recently we realized that it is 
quite impossible to describe soliton-like solutions of the following generalization of the Burgers 
equation (GBE from now on): 

TUtt + UUx + Ut- KUxx + JU{U^ - Zq) = 0, (17) 

using the ansatz ((H)), which is more general than the generalized Hirota-like ansatzes employed 
in [SJIH]. The first thing we are going to do in the situation when the ansatz-based methods fail, 
is to perform a qualitative study of the ODE system describing the TW solutions of the GBE 
and answer the question whether (and when) soliton-like solutions exist. 

Below we analyze in detail the whole set of TW solutions of equation (|17j) . The results of this 
analysis are used in the following section in constructing the approximated soliton-like solutions. 
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Inserting the ansatz u = U{£,) = U{x + fit) into the equations (fT7|l we obtain the second order 
ODE 

h^ + ^(^' + U)+lU {U' - zl) = 0. (18) 
h = r/i^ — K. The scaling transformation 



zq \ n 

leads us to the equation 

-l2 



^ + A^{f, + U) + U{U'-l)=0, (19) 



where A = 1/ \/hrf^ and fi = h/zq (we omit bars over the variables in the following formulae) . 
Let us rewrite H19|) as the following system of the first order ODEs: 

^ = -W, ^ = U{U'-l)-AW{f. + U). (20) 

System (|2U|) has three stationary points: -Bo(0, 0), Bi = (—1,0) and B2 = (1,0). Analysis 
shows [T3] that the critical point Bq is always a saddle. Depending on the values of the para- 
meters, two other critical points are either nodes or foci. As was declared earlier, we are looking 
for the homoclinic trajectory and since the dynamical system (|20() is of dissipative type, this 
trajectory can exist only for special values of the parameters. In our attempts to capture the 
homoclinic loop we begin with studying the condition of the limit cycle creation in proximity of 
the critical point i?2- If there exists a limit cycle lying in the right half plane and surrounding 
the critical point B2, then its evolution in the vicinity of the unmovable saddle point Bq, caused 
by the change of the driven parameter n, most probably will lead to the homoclinic bifurcation. 

It can be shown by inspection of the Andronov-Hopf theorem conditions and their fulfillment 
jl4l I13j that in vicinity of the critical value of the parameter fi = —1 an unstable limit cycle is 
created with zero radius. A further evolution of system (|2()|1 was undertaken by means of the 
numerical simulation in which we put A = 1. It shows that the radius of the periodic trajectory 
grows with the parameter fi, until the homoclinic bifurcation takes place at /i ~ —0.836. Note 
that the homoclinic solution intersects the horizontal axis at the point [/ = ~ 1.426095. 
On further growth of the parameter fi the homoclinic loop is destroyed and one of the phase 
trajectories starting at the unstable saddle Bi reaches the critical point -82- Patterns of the phase 
trajectories of system (|20|) obtained by means of numerical simulations are shown in Fig. |2j 

Thus, the homoclinic loop exists in system (|20|) merely at the specific values of the parameters. 
The construction of the approximated soliton-like solution, undertaken below is based in essential 
way upon the preliminary information obtained in this section. 

3.2 Approximation of the soliton-like TW solution to the GBE 

The results of qualitative analysis tell us that equation (|17|) possesses the SW solution for the 
specified values of the parameters. Let us assume that it can be described by means of the 
formula 0- Using the fact that the SW solution is represented by the homoclinic trajectory 
of system (fT^. that is bi-asymptotic to the origin, we can use for its description the following 
representation 

m 

J2 ak exp [ka^] 
1+ ^ br exp [ra^] 

r=l 
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that gives us proper asymptotics when ^ — > itcxD. The problem is that the formula (|21|1 does not 
describe the SW solution whenever m is a finite number, as it is stated below. 

Lemma 1. A solitary wave solution of the equation ()17() cannot be described by the formula (|21|). 

Proof. Since there is a one-to-one correspondence between the solutions of equations (|18j) 
and ()19() . we consider the latter one, trying to represent its homoclinic solution by means of the 
formula 

m 

Ofc exp [kaT] 

^(^) = • (22) 

1 + ^ 6^ exp [raT] 

r=l 

Let us denote the numerator and denominator of the above expression by F{T), G{T) respec- 
tively. Inserting the ansatz (|22|l into the equation (fT9|) and multiplying the resulting equation 
by , we get the expression 

G{F"G - G"F) + [A{F + fiG) - 2G']{F'G - G'F) + F{F^ - G^) = 0. 

Now, gathering the coefficients standing at functions exp (aT), exp [(3m + 2)aT] and assuming 
that the coefficients ai, am, b^+i are nonzero, we obtain the following system of algebraic 
equations: 

+ Afia — 1 = 0, — Afia — 1 = 0, 

that leads to a contradiction when A is nonzero. ■ 
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Of course, the formula ()22() does not give the most general expression for a solution, having 
the proper asymptotics at T — > ztco and belonging to the family (jSJ. For example, we can take 
the formula 

m+ni 

ak exp [kaT] 

U{T) = ^^^i— (23) 

1 + Yli exp [raT] 

r=n2 

instead of (|22|) for any natural numbers ni, ?i2, n^, and m. But our experiments with this formula 
in the cases when the numerator and denominator contained, respectively, two and three terms, 
gave a result which is identical to that stated in Lemma^ We do not mention it here because it 
is too cumbersome. Concerning an expression containing more terms the result is unclear as yet, 
but it is rather doubtful that we would obtain a positive answer by incorporating extra terms, 
for a number of extra algebraic equations accompanying the "inflation" of the formula H23|) 
exceeds in an essential way a number of extra coefficients a^, hr to be determined. 

So putting aside attempts to discover the analytical expression for soliton-like solutions of 
the GBE, we are searching for approximated solutions to the equation (|19j) . represented by the 
series (dH), 



U{T) = I 



an exp [aTn] when T < 0, 



(24) 



oo 



hn exp [—(3Tn] when T > 0, 



V n=l 

with the additional condition 

[/(O) = x„ C/'(0) = 0, (25) 

where x^, is a point of intersection of the homoclinic trajectory with the horizontal axis, which, 
generally speaking, can be obtained from the numerical experiment. With this conditions the 
two branches of the homoclinic trajectory are "tied up" in the point (x*, 0). 

There are two ways of determining the unknown constants and a. Firstly, it is possible to 
insert the series 

oo 

^{T) =Y,ake^v[aTk\ (26) 

k=l 

into (|17j) and equate to zero the coefficients standing at different powers of the function exp [aT]. 
This way we obtain the recurrence a^+i = ''Pk+i{(i-i^^2^ ■ ■ ■ jOfc)- Another possibility is to use 
the initial conditions (|25j) and their differential consequences. As in the previous case, this 
approach is quite algorithmic for any second order ODE that is linear with respect to the highest 
derivatives. The approach is based on the following statement, which can be easily proved by 
means of mathematical induction. 

Lemma 2. Consider the Cauchy problem 

U"{0=F{^,U,U'), (27) 
U{0) = xo, U'{0) = xi. (28) 

Suppose we are looking for the approximate solution to ((77|l . (P5|l in the form 

N 

U{0 = ^akexp[a^k]. (29) 

k=l 
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Then the coefficients {ak}j.^i can he obtained from the conditions and the N — 2 equations 
N 

fe=l 

J = 0, 1, . . . , A^ — 3, where denotes the total derivative. The functions <l>j(C/(0), . . . , U^^^^\Q)) 
are recurrently expressed merely by the initial conditions xq, xi and do not depend upon a^- 

So according to LemmaI21 using Cauchy data (PH)) and their differential consequences we get 
the hnear system 



(« • 1)° 


(« • 2)° 


(a.3)0 •• 


• (a • Nf ' 




fll 




Xo 


(a -1)1 


(a -2)1 


(a -3)1 •• 


{a ■ N f 








Xi 


(a -1)2 


(a • 2)2 


(a -3)2 • • 


■ (a • Nf 




03 




$0 


(a -1)^-1 


(a • 2)^-1 


(a -3)^-1 •• 


■ {a-N)^~\ 




_ ON _ 







The matrix Mjy standing at the RHS of equation (|30j) is of the Vandermonde type ^Hl- Its 
determinant is given by the formula 

JV(JV-l) 

detMAT = a^~ • 0! • l! • 2! • • • (A^ - 1)!, 
so it is always nonsingular and the system has a unique solution for each A^. 

Remark. The above lemma says nothing about the convergence of the series (|29|) to the real 
solution of the initial value problem ()27p. (|28() as A^ — > oo. This problem needs a special 
investigation. 

Our analysis suggests that the upper and lower branches of the homoclinic trajectory are 
so strongly asymmetric that the convergency of the approximated solution based on the first 
method of determining the coefficients of series expansions (|24() takes place merely for the lower 
one, whereas the upper branch should be approximated with the second method only. 

Thus, we approximate the lower branch of homoclinic solution by the series 

oo 

V'(r) = J^6„exp [-/3rn], 

n=l 

truncated on a number A^. The finite series is inserted into the factorized equation (|19() and 
the coefficients standing at the subsequent powers of function exp[— /3T] are equated to zero. 
For those values of the parameters that have been employed in the Section 2, the coefficient 
at {exp[— /3T]}^ nullifies if /? « 0.6658. Equating to zero the coefficients of {exp[— /?r]}^ , we 
obtain the recurrence for the coefficients 6^, k = 2,3, . . . , N , while bi remains unidentified. 
We omit the presentation of this recurrence since it is irregular and cannot be expressed it 
in the analytical form. Therefore all the calculations were performed with the help of the 
"Mathematica" package. They were based in an essential way upon the assumption that the 
phase trajectory approaches the point U = x^ = 1.426095 of the horizontal axis as T — > 0. 
Approximate solutions corresponding to different values of the number A^ are shown in Fig. |31 
Unfortunately we fail to employ a similar series decomposition and the method of finding the 
unknown coefficients for the upper branch, since the series do not converge to the numerical 
solution in proximity of the far end of the homoclinic trajectory, as it is seen in Fig. 0J In fact, 
the series is convergent yet the limiting function does not coincide with the solution, obtained by 
means of the highly precise numerical scheme. So in order to determine the series coefficients. 
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Figure 3. The lower branch of the homoclinic 
loop of system (jJSJ: numerical solution (solid) 
and the truncated series (|24|) (dashed). The 
coefficients Ofc are defined by means of the first 
method. The bold dashed line is obtained for 
N = 40. 




Figure 5. The upper branch of the homoclinic 
loop of system (|20|) : numerical solution (solid) 
and the truncated series (dashed). The 

coefficients are defined by means of the sec- 
ond method. The bold dashed line corresponds 
to TV = 20. 




0.2 0.4 0.6 0.8 1 1.2 1.4 

Figure 4. The upper branch of the homoclinic 
loop of system (|2()|l : numerical solution (solid) 



and the truncated series (|24|l (dashed). The 
coefficients are defined by means of the first 
method. 




-4 -2 2 4 6 

Figure 6. Solitary wave solution of equation (|17|l 
(solid) and its approximation (dashed). 



we turn to the second scheme. Keeping in the series H26() elements and substituting the 
truncated series into equation (|19)) . we obtain from (|25j) two conditions for the coefficients ak- 

N N 

U{0) = "^ak = x^ = 1.426095, U'{0) = ^ /ca^ = 0. (31) 

k=l k=l 

The first differential consequence of these conditions gives us equation (|19() itself. Taking it into 
the point T = and employing equations (|31() . we obtain: 

N 

U"{0) = Y^iakfau = (1 - xl)x. = h{x,). 

k=l 

In order to define all the coefficients of the series truncated on A^-th element we need — 3 
extra conditions. We can obtain these missing conditions from further differential consequences, 
as it is stated in Lemma |21 Coefficients {ak\k=i satisfy the linear system PU)) that delivers 
a unique solution, providing the RHS and the parameter a are determined. The remaining 
N — 3 components of the vector, positioned at the RHS of system PU]) . were obtained by 
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means of some recurrent procedure, written in "Mathematica" package. In order to obtain the 
parameter a, we insert the solution (|26|) into (|19|) and equate to zero the coefficient a'^ +Afj,a — 1, 
at the lowest power of the exponential function. For A = 1 and /U = —0.836 this gives us the 
value a = 1.5018. 

The above procedure was used to obtain the upper branch of the homoclinic curve of sys- 
tem H20() (or what is the same, the left "wing" of the SW solution of the equation ()19() ). The 
results obtained are shown in Fig. [3 The upper and lower branches sewed together are shown in 
Fig. 121 from which we can see that suggested approach quite well approximates the homoclinic 
solution of equation (fT7)). 

4 Concluding remarks 

We have reviewed in this study the methods and tools that allow one to describe the solutions 
of the nonlinear PDEs with the required features and asymptotic behavior. The exact solutions 
presented here are obtained with the help of the ansatz which proves to be effective when 
the methods described in 13IH] fail. We concentrated mainly upon the soliton-like TW solutions 
and their description, yet similar tools can be applied to the description of kink-like, periodic 
and multi-periodic solutions as well. 

In contrast to the classical self-similarity methods, ansatz-based methods (3 IH1 El ^1 enable 
to obtain the exact solution with required properties an asymptotics (such e.g. as the periodic, 
soliton-like and kink-like TW solutions). This is due to the fact that within these methods 
solutions are defined as algebraic combinations of certain set of functions (elementary of special 
ones), closed with respect to differentiation and algebraic operations, being represented in a form, 
enabling to control some features and the asymptotic behavior of the solutions to be described. 
Yet efficient employment of the direct algebraic balance methods is related to one's ability to 
solve large systems of nonlinear algebraic equations. In addition to this purely technical obstacle, 
there is a number of evidence that none of the already known version of the balance algebraic 
method is fully universal. For this reason the authors put forward the method of the asymptotic 
description of the soliton-like TW solutions by an infinite series of exponential functions. The 
effectiveness of this method, connected with some general properties of the set of exponential 
functions |T71, is demonstrated by its application to the GBE. Previous qualitative investigations 
showed the existence of such solutions for the specific values of the parameters and further on it 
was shown that these solutions cannot be described as finite algebraic combination of exponential 
functions. 

Let us say a few words about the problem of the truncated series (|24|) convergence to the 
exact solution. A uniform convergency of a series of exponential functions to an exact solu- 
tion can be easily stated e.g. for the nonlinear wave equation, because the expressions for the 
coefficients a„ and 6„ are known in this case J3]- Unfortunately, for GBE the recurrent expres- 
sions for indices a„, 6„ of the decomposition H24|l are not known and the problem of the series 
convergence remains still open. 

The method of approximation of solitary wave regimes, described in Section 3, is relatively 
easy to perform in case when the factorized system is two-dimensional. In the less trivial cases, 
when the factorized system is multidimensional, the homoclinic trajectories representing the 
soliton-like solutions become much more complicated. As an example let us mention a Shil'nikov 
homoclinic trajectory associated with the saddle-focus of a three-dimensional dynamical sys- 
tem |18| IH]. In order to deliver an approximated description of such trajectory following the 
methodology outlined in Section 3 of this study, it is necessary to incorporate wider class of 



During the preparation of this proof, we wre pointed to the paper |16j in which similar ansatzes were apphed 
to a family of weakly nonlinear equations of Boltzmann type. 
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functions. In earlier attempts at describing the Shil'nikov homoclinics a combination of 
the exponential and trigonometric functions was used. The approximation employed in |19j was 
not very effective because of the great technical difficulties encountered. Since then a num- 
ber of powerful programs for the symbolic calculations were developed, and it is our hope that 
approximation similar to that used in actually can be performed technically much more 
effectively. 
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